LP-INTEGRABILITY OF THE GRADIENT OF SOLUTIONS TO 
QUASILINEAR SYSTEMS WITH DISCONTINUOUS 

COEFFICIENTS 
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Abstract. The Dirichlet problem for a class of quasilinear elliptic systems of 
equations with small-BMO coefficients in Reifenberg-flat domain is considered. 
The lower order terms supposed to satisfy controlled growth conditions in u and 
Du. ft is obtained LP-integrability with p > 2 of Du where p depends explicitly on 
the data. An analogous result is obtained also for the Cauchy-Dirichlet problem 
for quasilinear parabolic systems. 

1. Introduction 

In the present work we study the integrabihty properties of the weak solutions of 
the following Dirichlet problem 

D aiji'^-' {x)D pu^ {x) + af{x, u)^ = hi{x, u, Du) a. a. x eVL 
vl{x) =0 on dVt 

where Vt C M", n > 2 is a bounded Reifenberg flat domain (see Definition 2). 
The principal coefficients are discontinuous with "small" discontinuity expressed 
in terms of their bounded mean oscillation (BMO) in Q (cf. [20]). The matrix 
A(x) = {Aif{x)}if^j;^ verifies 

(2) > ^1^1' V^GM^x", ||A|U,f. < M 

with some positive constants A and M. The non linear terms 

a(x, u) = {a"(x, u)}"<^^ and b(x, u, z) = u, z)}j<Ar 

supposed to be Caratheodory functions for x G fi, u G M^, z G M^^" and satisfy 
controlled growth conditions. Namely, for |u|, |z| — t- oo we have 



(1) 



(3) a"(a;, u) = 0((y9i(a;) + |u|"+2) and bi{x,u,z) = 0{ip2{x) + \u 
with (^1 G LP{Q), p > 2 and (^2 G L'^i^), Q > 



n + 2 , 
n-2 
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Our aim is to show tliat tlie problem (1) satisfies the Calderon-Zygmund prop- 
erty when Q is (5, i?)-Reifenberg flat and the coefficients are {S, i?)-vanishing in Q. 
Precisely, each bounded weak solution u e VFq'^ fl L°°(n;R^) of (1) gains better 
regularity from the data (pi and (p2 and belongs to Wq'"^^^^'^ ^ fl L°°{fl;R^) where 
q* is the Sobolev conjugate of q (see (4)). 

Similar result is obtained also for the Cauchy-Dirichlet problem for the parabolic 
quasilinear system 

iul - Da(Aif(x, t)Dpu^ + af (x, t, u)) = hi{x, t, u, Du)) a.a. (x, t) e Q 
[u{x,t)^0 {x,t)edQ 

in a cylinder Q = fl x {Q,T) where fl is (5, i?)-Reifenberg flat and dQ = flL) {dfl x 
(0,T)} is the parabolic boundary. 

The problem of integrability and regularity of the solutions of linear and quasilin- 
ear elliptic/parabolic equations and systems is widely studied. Let us start with the 
classical results concerning equations/systems with smooth coefficients presented in 
the monographs [24, 25]. In the scalar case, N = 1, the notorious results of De 
Giorgi [10] and Nash [30] assert Holder regularity of the solutions of linear diver- 
gence form equations with only L°° principal coefficients. One remarkable result 
that permits to obtain higher integrability of the weak solutions is due to Gehring 
[16]. He studied integrability properties of functions satisfying the reverse Holder 
inequality. It was noticed that some power of the gradient of the weak solutions 
satisfies local reverse Holder inequality. Modifying the Gehring lemma, Giaquinta 
and Modica [18] firstly obtain higher integrability of solutions of divergence form 
quasilinear elliptic equations. For the sake of completeness we give this result as it 
is presented in the monograph by Giaquinta [17, Theorem V.2.3]. 

Theorem 1. Suppose that g e Li{Q), F e L«+'^(Q), g, F > 0, q > 1, 5 > and 
■I g'^dx < B (-f gdx] + / F^dx + 6-1 g'^dx 

JBr{x) \JB2r(x) ) JB2r{x) JB2r{x) 

for a.a. x e Vl, R < ^ min{(i(a;, 90), i?o} where Rq > 0, B > 1, 9 e [0,1). Then 
g e LP'^^'^in) and 




for any ball B2R C fl, 2R < Rq, p G [q,Po) where C > 0, po > q depend only on 
B,e,q,n. 
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There are various generalizations of the above theorem permitting to study elhptic 
and parabohc problems with Dirichlet and Neumann boundary conditions (see [2, 
3, 4, 5, 11, 15, 27]). Other results concerning higher integr ability of divergence 
form quasilinear equations and variational equations could be found in [9, 17, 19, 
21, 28]. The L^-estimates of derivatives obtained such way laid the foundation to 
the so-called "direct method" of proving partial regularity of solutions. Recently, 
the method of A-harmonic approximation permits to study the regularity of the 
solutions without the use of the Gehring lemma. For more details we refer the 
reader to [1, 12, 13, 14], see also the references therein. 

The regularity theory for linear operators with smooth data was extended on 
operators with discontinuous coefficients defined in rough domains. In [6, 7, 8] the 
authors consider divergence form elliptic and parabolic equations and systems with 
BMO coefficients in Reifenberg flat domain with Dirichlet boundary conditions 
extending such way the known results on operators with VMO coefficients too (see 
also [28, 34, 35, 22, 23] and the references therein). In [15] a reverse H51der inequality 
is established for quasilinear elliptic systems with principal coefficient being VMO in 
X and under controlled growth conditions over the lower order terms. It permits the 
authors to obtain interior Holder continuity of solutions to scalar equations as well as 
partial Holder regularity of solutions to systems. In [31, 32] global Holder regularity 
of solutions to elliptic quasilinear equations with VMO in x principal coefficients 
is proved under strictly controlled growth conditions. Later this result is extended 
for quasilinear elliptic and parabolic equations in Reifenberg flat domains supposing 
controlled growth conditions and Dirichlet boundary data (see [11, 33, 36, 37]). 

In the present work wc extend the results from [37] to elliptic and parabolic 
systems with discontinuous data. Making use of the linear L^-theory for systems, 
developed in [7, 8] and the bootstrap method we prove Du G with r depending 
explicitly on the data ipi and ip2 in (3). 



2. Elliptic systems, definitions and main result 

In the following we use the standard notations: 

• X = (xi, . . . , a;„) G M", p > and Bp{x) = {y e W : \x - y\ < p}. 

• let Q C IR" be a bounded domain, x E Q and denote ilp{x) = Q.r\Bp{x). 

• M^^" is the set oi N x n-matrices. 

• For a vector function u — {u^,. . . , u^) : — >■ IR-^ we write 



|u|2= ^ |«^f, D^u^^ 

j<N 



_d_ 

dxn 



-u 
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a<n 

j<N 

• Let / : r2 ^ M and be the Lebesgue measure of ^2, then 

• For u e LP{n; R^) write || u\\p^Q instead of ||u||ip(!^.]Rjv). 

• For each s e (1, oo) recall that s* means the Sobolev conjugate of s 

(4) .* = /^ 

I arbitrary large number > 1 if s > n. 

For the function spaces we follow the notions of the monographs [24, 28]. Through all 
the paper the standard summation convention on repeated upper and lower indexes 
is adopted. The letter C is used for various constants and may change from one 
occurrence to another. 

In [38] Reifenberg introduced a class of domains with rough boundary that can 
be approximated by hyperplanes at every point and at every scale. Namely 

Definition 2. The domain Q is (5, i?)-Reifenberg flat if there exist positive constants 
R, 6 < 1 such that for each x E dQ and each p e (0, R) there is a local coordinate 
system {yi, . . . , y„} with the property 

(5) Bp{x) n {yn > 5p} C Qp{x) C Bp{x) (1 {y^ > -5p}. 

Reifenberg arrived at that concept of flatness in his studies on Plateau's problem 
in higher dimensions and he proved that such a domain is locally a topological disc 
when 6 is small enough. It is easy to see that a C^-domain is a Reifenberg fiat with 
5 — > as it! — > 0. A domain with Lipschitz boundary with a Lipschitz constant less 
than S also verifies the condition (5) if S is small enough (say S < 1/8). But the 
class of Reifenberg's domains is much more wider and contain domains with fractal 
boundaries. For instance, consider a self-similar snowflake Sj^. It is a flat version of 
the Koch snowflake Sj^/^ where the angle of the spike with respect to the horizontal is 
13. A domain f2 C with Sf^ C dVl is a Reifenberg fiat if < sin /3 < 5 < 1/8. This 
kind of flatness exhibits minimal geometrical conditions necessary for some natural 
properties in analysis and potential theory to hold. For more detailed overview of 
the properties of these domains we refer the reader to the papers [29, 39]. 

From (5) it follows that dQ satisfies the (A)-property (cf. [9, 17, 24]). Precisely, 
the measure \^p{x) \ is 5-comparable to \Bp{x) \ , that is there exists a positive constant 
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A{6) < 1/2 such that 

(A) A{6)\B,{x)\ < \Q,{x)\ < (1 - A{6))\B,{x)\ 

for any fixed x G dQ, p G (0, i?) and S G (0,1). This condition excludes that Q 
may have sharp outward and inward cusps. Moreover, for small S they can be 
approximated in a uniform way by Lipschitz domains with a Lipschitz constant less 
then S (see [8, Lemma 5.1]). As consequence, they are '^-extension domains, 
1 < p < oo, hence the usual extension theorems, the Sobolev and Sobolev-Poincare 
inequalities are valid in Q,. 

To describe the discontinuity of the principal coefficients we need of the following 

Definition 3. We say that a function a{x) is a {6, R) -vanishing if there exist positive 
constants R and S < 1 such that 

(6) sup sup/ \a{y)-aa/^)\^dy<5^, aa{x) = f a{y)dy. 

We suppose that all A'^^^{x) are (5, i?)- vanishing. It implies that A G BMO{Q) 
with a small BMO norm || A||* < S. 

The nonlinear terms a(x, u) and b(x, u, z) are Caratheodory functions for x G il, 
u G ]R^, z G M^^" and satisfy the controlled growth conditions 

(7) |a(a;,u)| < A((^i(x) + |u|^), cpi e LP{Q), p > 2 

2n 
n + 2 

with some positive constant A. In the particular case n ~2 the powers of |u| could 
be arbitrary positive numbers while the growth of |z| is quadratic (cf. [17, 24]). 

Under a weak solution to the problem (1) we mean a function u G W^o'^(f^; M^), 
1 < p < oo satisfying 



(8) |b(a;,u,z)| < A((^2(a;) + |u|^ + Iz]"^), ip2 e L'^iVt), q > 



f Af{x)D^u^{x)D^x\x)dx+Y. [ at{xM^))DaX\x)dx 



a<n ' 

i<N 



+ ^ bi(x, u{x), Du{x))x\x)dx = 0, j^l,... ,N 

for all X ^ W^o^'^ (fi; R^), p' = p/{p — 1). The conditions (7) and (8) are the natural 
ones that ensure convergence of the integrals above. Moreover, they are optimal 
since a growth of the gradient greater than leads to unbounded solutions as it 
is seen from the following example (cf. [25, 32]). The function u{x) G W^''^{Bi{0)), 

r — 2 

u{x) = \x\~ is a solution of the equation Au — C\Du\^ in Bi{0). Note that u{x) ^ 
L~(Bi(0)) for ^ < r < 2. 



6 



L.G. SOFTOVA 



In generally we cannot expect boundedness of each solution of (1) unless wc add 
some structural conditions. Consider, for instance, the system 

DaiA'^ix,u,Du)) = bi{x,u,Du) xen 

where 

l3<nj<N 

are measurable in x E Q. Assume a pointwise coercive and sign conditions, both 
of them for large values of the corresponding component of u, precisely: for every 
i G {1, . . . , N} there exist constants 6'*, M', u e (0, +oo) such that for > 6^ we 
have 

(9) < a<n 

< bi(x, u, for a.a. x e Q, e M^^". 

2n 

Suppose (7), (8) and (9) and let u e W^'^ n L^(n;M^) be a weak solution of (1) 
then for each i e {1, . . . , N} 

supu' <e' + K' 

n 

where depend on M\n. \ Q\ and u (sec [26]). 

Theorem 4. Let u G W^''^ r]L°°{Q; M^) be a weak solution of the problem (1) under 
the conditions (2), (7) and (8). Then there exists a small number 5q > such that if 
fl is {6, R)-Reifenberg fiat domain and Af^^{x) are {6, R) -vanishing with S < Sq < 1 

then 

(10) u G Wo'" n L^{n; M^) with r = mm{p, q*}. 

Proof. In [17, Chapter 5] Giaquinta considers quasilinear strongly elliptic systems 
with L°° principal coefficients, under the conditions (7) and (8). Making use of 
the reverse Holder's inequality and the version of the Gehring lemma it is shown 
that there exists an exponent rg > 2 such that u G Wi^^°{^l,'M.^) (cf. [17, The- 
orem V.2.3],[9, Chapter III] or [28, Lemma 3.2.23] ). Since, roughly speaking, 
Caccioppoli-typc inequalities hold up to the boundary, the method for obtaining 
higher integrability can be carried over up to the boundary. In [17, Chapter 5] it 
is done for the Dirichlet problem in Lipschitz domain. Since the Reifenberg flat 
domain can be uniformly approximated by Lipschitz domains the same result still 
holds true. Precisely, there is ro > 2 such that 

(11) ||^u||,,f^<Ar VrG[2,ro) 
where N and ro depend on n. A, A, ||</?i||p,f2, ||</?2||g,n, |^^|, ||-Du||2,f2. 
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Let n > 2 and u G W]'"^ n L°°(1];M^) be a solution of (1). Fixing that solution 
in the nonlinear terms we get the linearized problem 



Da{A"f {x)Di3U^) — fi{x) — div(Aj(x)) a.a. x E fl 
u(x) = on do, 



(12) 
where 

fi{x) =bi{x, u, Du), f (x) = h{x, u, Du), 
Ai{x) ={al{x,u), . . . ,a"(a;, u)), A(a;) = (Ai(a;), . . . ,An{x)) 

and by (7), (8) and (11) we get 

\\M\p,n<C^y4,,n+\\u\\^,^ 

( n + 2 
ll<^2||,i,n + ||u||;,-^+||L'u| 



(13) 



71 + 2 n + 2 

I qi(.n+2) 



with p > 2 and qi = min Further, for aU fi G L'^^{^), i = I,... ,N 

there exists a vector field ¥i{x) G L'^i(f2,M") such that fi{x) = divFj(a;). Denote 
¥{x) = (Fi(a;), . . . ,Fjv(a:)), then by [32, Lemma 3.1]) we have 

(14) mk,^<CmU,a, gi = min{g,-^}. 

Thus the problem (12) becomes 



(15) 



Da{Atf{x)Di3U^{x)) = div(Fi(x) - Ai{x)) a.a. ,t G Q 
u(x) = 0, X G Sil. 



For linear systems as above we dispose with the regularity result of Byun and Wang 

[8, Theorem 1.7] that asserts there exists a small positive constant S = S{X,p,n, N) 
such that for each (5, -R)-vanishing A'^^^ , for each (5, i?)-Reifenberg flat Q, and for 
each matrix function F — A G L''^ ((]; M^^"), with ri = min{p,ql}, the solution 
u G Wq'^ n L~(0;R^) of (15) belongs to PFq'"^ n L~(f];M^) and the following 
estimate holds 

(16) \\Du\\r,,n < C\\¥ - A\\r,,n, n = min{p, g*} 

with C = C(A,p,n,A^, \Q\). 

Our goal is to show the inclusion Du G U{0,;M^^'^) with r = min{p,q*}. For 
this we study the following cases: 

1) li q < ^ then qi — q in (14) and ri = r — min{p, q*}. 
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2) Ifg>^,thengi = ^and 



Ton Ton 

II < n 



Qi 



n + 2 — To n + 2 

rpn 

n + 2 



arbitrary large number > 1 if > n. 



Consider again two sub-cases: 

2a) If n > then ri = min{p, ^ ' If — P then the theorem holds 

true otherwise Du G UUn : M"^^) with ri = 7"°" . 

2{,) If n < then gj' is arbitrary large number, that implies ri = p and 
the theorem holds true once again. 

It is easy to see that ri = r unless 

ran ron 

< q and < n 



n+2 n+2 

when u e Wq'"'' D L°°(Q;K^) with ri = It holds for any solution of the 

linearized problem (15) including the one fixed in the coefficients in (1). 

Consider once again (13) with Du e L^i(Q;M^^"). Hence fi e L^^in) with 
q2 = min{g, („_)_2)2^"p(„_)_2) } and the associated vector- field Fj belongs to L'^'^{Vt] MJ^). 
Than Fj - e U^{Vt;W) with r2 = min{p,g^}. Applying [8, Theorem 1.7] to 
system (15) and repeating the same procedure as above we get that the theorem 
holds with r2 = r if 

(17) 

> n or 11] — -, — > q. 



' (n + 2)2 - roin + 2)- ' (n + 2)2 - ro(n + 2) 

Otherwise r2 = , , ' "7" , o\ if 

/ION '^0^^ J ^0^^ 



, , ^ < q and -r ; < n 

' ' (n + 2)2 -ro(n + 2) (n + 2)2 - ro(n + 2) 

Repeating the same procedure /c-times we get that the assertion holds if 



Ton 

(n + 2)*= - ro eS ^'(^ + 2)*^-' 
Direct calculations give that (19) is equivalent to 



k > min 



ro /, n + 2n r /•(,(2f7 + fy// + 2) / n + 2 



+ 1 



ro — 2' n -''L g(n + 2)(ro — 2) ' n 

where [x] means the integer part of x. 

The case n = 2 is simpler and is left to the reader. □ 
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3. QUASILINEAR PARABOLIC SYSTEMS 

Let Q = O X (0,T) be a cylinder in M'*+^ with Q. being (5, i?)-Reifenberg flat. 
Denote by Cp the parabolic cylinder 

Cp{x,t) ^Bp{x) X {t- p\t), Qp{x,t) ^Qr\Cp{x,t) for(x,t)eg, 



-f a(y,T)dydT = --— / a(y,r)dydi 



Let 1 < r < oo and u : Q M^. 

1. The space W^'^{Q]Ml^) consists of all functions u e U{Q]'MJ^) having a finite 
norm 

1 1 1 1 1 1 1 1 II 1 1 T* 

2. The space Wy{Q;R^) = L''{0,T;W^^''{n;R^)) f] W^'''{0,T;W-^''' {n;R^)), 
r' = r/(r — 1) consists of the functions u G W^'^{Q; M^) for which there exist vector 
functions g e L'-{Q;R^) and F e L''(g;M^^") such that 

Ut — divF — g a.e. in Q 

in the sense of distributions, that is, for each vector function x ^ C^{Q) with 
Xix, T) = holds 

(20) / u-xtdxdt^ I (F-Dx + g-x)dxdt. 

The space W^'^{Q;'R^) is endowed by the norm 

Miw^iQ) = Mw^'^^iQ) + I (/^ iFr + igr) ^ I 

where the infimum is taken over all F and g satisfying (20). The closure of C^{Q) 

o 

with respect to this norm is denoted by Wl''^{Q; R^). 

3. V2{Q; R^) stands for the Banach space of all functions u e W^'^iQ; R^) for which 

l|u|k2(Q;K^) = esssup ||u(-,i)||2,n + ||-Du||2,Q < oo . 



te[o,T] 

4. V2'^{Q]R^) consists of all u e V2{Q;R'^) that are continuous in t with respect 
to the norm of L'^{n;R^) 

hm ||u(.,i + At)-u(-,t)||2,n = 0. 

The norm in V2'^{Q', R^) is given by 

l|u||yM(Q.KiV) = max^ \\u{;t)\\2,n + \\Du\\2,q. 
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We consider the Cauchy-Dirichlet problem for the strongly parabolic quasilinear 
system 

(21) 

ul — Da{Aff DjjU^ + a^XjijU)) = bi{x,t,u, Du) a. a. {x,t) E Q 
u{x,t) = {x,t)edQ. 

The principal coefficients satisfy e L°°{Q) and 



(22) Al^(%>v\i\^ VeeM^><", i. = const>0. 



In addition we suppose that are (5, i?)-vanishing, that is 

(23) sup sup / \Af{y, r) - < Jdydr < 5' 

which implies small BMO norm of each A'^l^ . 

The functions af(a;,t, u), 6.4(2;, t,u,z) are Caratheodory ones and verify the con- 
trolled growth conditions (see [2, 25]) 



(24) K{^^tM\ < HM^,t) + |u|^) 

i<N 

(25) ^ \bi(x,t,u,z)\ < A(V'2(a:,^) + |u|'^ + |z|^)■ 
^<A?" 

with V'l e LP{Q),p > 2 and V'2 e Li{Q), q > 

o 

A vector function u e fl L°°(Q;M'^) is a weak solution to (21) if for any 
function x e W^P(Q; ^^), T) = we have 

/ u%x, t)xt{x, t)dxdt - I (A"/(x, t)Dfsu^{x, t) + a"(x, u))i:>aX'(2;, i)cixcii 

+ / 6i(a;, t, u, L'u)x*(a;, t)da;di = . 
JQ 



Theorem 5. Let u e Vl^^'^ n L°°(Q;R^) 6e a weak solution to (21) under the 
conditions (22) -(25). Then there exists a small positive constant 5q <1 such that if 
Q, is (5, R)-Reifenberg flat and A'^f are {5, R)-vanishing with Q < 5 < 5q then 

u e WY n L°°(g; M^) with r = min{p, q**} 

.Here ,"4^ 

I arbitrary large number > 1 if q > n + 2. 
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Proof. The higher integrabihty of the gradient of the solution follows by the modifi- 
cation of the Gehring lemma due to Arkhipova [4, Theorem 1] which is very efficient 
for the study of parabolic systems with controlled growth conditions in domains 
with boundary dD, satisfying a kind of (A)-property. Recently, similar result is ob- 
tained in [11] for domains having strongly Lipschitz boundary. Since the Reifenberg 
flat domain can be approximated uniformly with Lipschitz domains with a small 
Lipschitz constant [8, Lemma 5.1]) we have that there exists Tq > 2 such that for 

any solution u e Wl'^{Q;R^) holds 

\\Du\\r,Q<N Vre[2,ro) 
where Tq and N depend on the data of the problem and ||£)u||2,q. Take a solution 

u e Wl'^ n M^) of (21) and fix it in the lower order terms. Thus we get the 

linearized problem 

ul — Da{A°^f D/su^) = fi{x,t) -I- div Aj(a;, t) a.a. {x,t) G Q 
u{x, t) = on dQ 



(26) 
where 



Ai{x, t) = {al{x, t, u), . . . , <(a;, t, u)), A(,x, t) = {af{x, t, u{x, t))}g^ 
fi{x, t) = bi{x, t, u, Du), f ^ {fi{x,t),... , fN{x, t)). 
Making use of the conditions (24) and (25) we get 

EII^»IUq<c(II^i|Ip,q + I|u|L",qJ 

i<N 

, K+4 w+4 

gl(n+4) ] 

i<N "+2 

with p > 2, gi = min {g, ^^^g^}. Let ft' C W be C^-domain such that Q ft' and 
consider the cylinder Q' = Q' x (0,T). Suppose that fi{x,t) is extended as zero out 
of Q. It is well known (cf. [25]) that for each /j e L'^^{Q') the linear problem 

'Fi - AF' = Fi - DaiS-^D^F') = /,(x, t) a.a. (x, t) e Q' 
F'(x, t)^0 on dQ' 

o 

has a unique solution G W^'_^{Q') and the estimate holds 

Denote by F = {F\ . . . ,F^), = {Aff - S'^^Sij)D^F^ and = (F^, . . . ,F,'^). 
Since the Sobolev trace theorem holds for domain with Reifenberg flat boundary we 



(27) 
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get that F'^ls E W ''i' (S) where S = dQ x (0,T). Consider the hnear system 

\u'-F%-D^{A^fD0{a^ -F^)) 

= D^iiA^f - 6"f'6,j)Df^F^) + div (A,) 
= div {Wi + Ai) t)eQ 
u - F = -F on dQ. 

By the imbedding theorems, DF^ G L'l* (cf. [25, Ch. II, Lemma 3.3]) and 

Hence Fj + Aj e L^'^^iQ) with ri = mm{p,ql*}. Applying [7, Corollary 2.10] on the 
linear problem (27) we get 

\\D{n - F)\U,Q < C(l + + ||A|U) < C(l + ||f + ||A|U). 

Hence 

|PU||^,,Q < C(l + ||f + ||A||p,q). 

Applying the bootstrapping arguments , we obtain as in the elliptic case that the 
theorem holds after k iterations with 



k > min 



log(ro/(ro - 2)) ^ r^pg + 4) - 2(n + 2)] /^^^ n + 4 



log((n + 4)/(n + 2))J' L ° _ 2)(n + 4) / ° n + 2 



□ 



Acknowledgments. The author is indebted to the referee for the valuable remarks 
which have lead to improvement of the article. 
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